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iBSTBiCT 


The  conputation  methoJ  developed  in  this  thesis  proceeds 
from  the  theory  developed  by  Jones  [Ref.  1].  His  final 
equations  are  first  rewritten  in  matrix  format.  They  are 
then  organized  into  computaticnal  sequences  tliat  may  be 
translated  into  computer  programming  language  to  calculate 
the  aerodynamic  parameters  of  wings  of  arbitrary  planform. 


— '-i 


3 


TABLE  OF  CCNTENTS 


I.  BACKCFOUND . 7 

II.  INTRODUCTION . 10 

III.  MATRIX  NOTATION . 11 

A.  CIRCULATION  AND  PRESSURE  RELATIONS  .  11 

B.  CALCULATION  OF  DERIVATIVES  .  13 

1.  Chorduise  Derivatives . 13 

2.  Spanwise  Derivatives  . . 17 

C.  INTEGRAL  ECUATION  FOE  HING  SLOPE . 20 

1.  First  Wing  Slope  Function . 20 

2.  Second  Wing  Slope  Function . 23 

3.  Final  King  Slope  Function . 24 

IV.  SUMMARY  OF  KEY  RELATIONSHIPS . 31 

A.  PRELIMINARY  CALCULATIONS  .  31 

B.  CALCULATION  OF  THE  A  AND  B  MATRICES . 32 

1.  The  A'  and  B'  Matrices . 32 

2.  The  V  Matrix . 34 

3.  The  A  and  3  Matrices . 35 

C.  AritITIONAL  LIFT . 35 

D.  BASIC  LIFT . 36 

V.  CCNCLUniONS . 38 

LIST  OF  FEPEEENCES . 3  9 

BIBLIOGRAPHY  .  40 

INITIAL  DISTRIBUTION  LIST  .  41 


4 


LIST  OF  TABLES 

1.  Elfitent s  of  Da*  and  Db*  Along  the  Three 


Prii -ipal  Diagonals  .  27 

2.  Elements  of  E*  Along  the  Four  Principal 

Diagonals . 29 


5 


LIST  OF  FIGURES 


3.1  Index  notation . 12 

1.2  Matrix  Da . 14 

’;.3  MatiixPb . . . 16 

!.4  Matrix  r: . 17 

3.5  2  ircu  idt  i.  on  point  . . 19 

3.6  Matrix  V . 24 

1.7  Dia^icnal  structure  of  Da'^and  Db** . 26 

1.8  Did'jonal  Structure  of  E* . 28 


6 


I 


I.  BACKGEOOND 

The  cvcrail  research  of  which  this  thesis  forms  a  part 
has  a  three-fold  purpose.  Firstly,  it  is  intended  to 
instruct  tl.e  student  about  the  f undement al  equations  that 
qovern  the  aerodynamics  of  wings  in  incompressible  flow. 
Secondly,  it  teaches  the  student  how  to  convert  the  funde- 
mental  equations  into  a  form  suitable  for  numerical  calcula¬ 
tion  on  the  digital  computer  and  to  carry  out  the  complex 
programming  required  for  this  purpose.  Thirdly,  once  the 
first  two  objectives  have  been  accomplished,  the  computer 
program  itself  provides  a  useful  pedagogical  and  design  tool 
for  illustrating  the  effects  of  various  wing  design  parame¬ 
ters  cn  the  final  aerodynamic  performance  of  a  wing. 

Lcdr.  J.  L.  Parks  made  a  creditable  start  toward  the 
first  objective  in  his  thesis  [Ref.  2J.  However,  in  the 
limited  time  available  to  him,  he  was  unable  to  make  signif¬ 
icant  progress  toward  the  second  or  third  objectives. 

It  was  evident  early  that  Parks*  computational  technique 
had  some  flaws.  The  expectation  was  that  these  would  be 
readily  discovered  and  soon  ccrrected  and  that  objectives 
two  and  three  would  be  rather  quickly  and  easily  completed. 
It  has  not  worked  out  that  way! 

As  the  investigation  actually  developed,  a  whole  scries 
of  obstacles  were  encountered  and  each  was  in  turn  eventu¬ 
ally  overcome.  But  the  process  was  neither  quick  nor  easy. 
In  every  case,  it  was  found  that  the  seeming  obstacle  was 
really  attributable  to  some  ccnceptual  error.  Once  the 
error  was  found,  the  obstacle  disappeared  and  the  investi¬ 
gator  gained  a  deeper  insight. 

Among  tJie  obstacles  that  were  encountered  and  eventually 
overcome  were  the  following; 
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1)  Confusion  about  the  proper  type  of  computational 
grid  to  use,  whether  staggered  or  unstaggered,  how 
fine,  and  so  on.  He  have  now  settled,  for  good 
reasons,  on  an  3  by  8  unstaggered  grid  over  the 
semi-span. 

2)  Confusion  about  how  to  deal  numerically  with  the 
singularity  that  occurs  in  the  governing  integral 
equations.  One  option  is  to  try  to  avoid  the 
problem  by  staggering  the  grid  of  the  field  points 
with  respect  to  the  grid  of  the  control  points. 
Instead  of  this,  we  have  resolved  the  indeterminacy 
by  rigorous  analysis  and  as  one  consequence,  are 
able  to  employ  a  simple  unstaggered  grid. 

3)  Confusion  about  how  to  evaluate  the  partial 
derivatives  of  the  circulation  function,  whether 
analytically  or  by  finite  difference  formulas.  It 
was  found  that  analytical  differentiation  is  incor¬ 
rect  and  that  finite  differences  must  be  used. 

4)  Confusion  about  the  validity  of  representing  the 
circulation  by  a  Fourier  series.  While  this  proce¬ 
dure  is  widely  advocated  in  the  technical  litera¬ 
ture,  we  found  that  a  much  simpler  and  clearer 
formulation  can  be  obtained  otherwise. 

5)  Confusion  about  the  boundary  conditions.  Special 
conditions  apply  to  the  leading  and  trailing  edges, 
the  wing  tip  and  at  midspan.  These  quite  complex 
boundary  conditions  have  now  been  fully  and  rigor¬ 
ously  analyzed  and  incorporated  consistently  into 
ouj  formula ticn  of  the  problem. 

If  ii;  evident  from  the  foregoing  discussion  that  this 
cor.  'M.ch  h  's  ’.counted  to  a  majcr  education  in  basic  aerody- 
naniL'  and  in  numerical  methods.  In  these  respects  it  has 
teen  a  deeply  rewarding  experience.  However,  in  view  of  the 
for  jjoing  td  ;:?  acles  and  problems,  the  time  schedule  has  of 
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course  been  greatly  delayed  frcm  that  which  was  initially 
anticipated.  Thus  there  are  nc  final  numerical  results  at 
this  particular  stage  of  the  investigation  and  there  is  no 
more  time  available  to  the  present  research  team.  Hence, 
this  final  aspect  will  have  to  be  completed  by  some  suLseg- 
uent  investigators.  Nevertheless,  what  the  present  effort 
has  produced  to  begueath  to  any  subsequent  worker  is  a 
sophisticated  and  refined  numerical  method.  This  method  has 
now  evolved  to  the  point  that  it  can  be  confidently  expected 
to  produce  a  reasonable  and  accurate  final  result. 
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II.  inteo^ction 

This  thesis  is  the  development  of  a  computational  method 
for  determining  the  aerodynamic  parmeters  of  wings  in  incom¬ 
pressible  flow,  from  the  theory  developed  by  Jones  [Bef.  1]. 
In  conjunction  with  Jones'  theory  this  analysis  has  solved 
the  problems  requiring  further  attention  pointed  out  by 
Parks  [Bef.  2].  Primarily  a  new  approach,  defining  the 
circulation  directly  instead  of  using  a  Fourier  series 
representation,  is  employed.  This  along  with  a  direct 
method  of  dealing  with  the  singularities  that  are  encoun¬ 
tered  in  ti;e  joverning  integral  equations  when  a  control 
poi  it  and  fi'^ld  point  coincide,  allows  the  computational 
jri ’  ■  fo:  th  v;e  points  to  be  identical.  This  analysis  is 
int  ‘.'ied  te  complement  that  of  Jones  by  taking  his  final 
theoretical  equations  and,  where  applicable,  presenting  them 
In  iiiitrix  t.  The  equations  are  then  organized  into  a 
compu tatioi al  sequence  which  can  be  further  translated  into 
programming  language.  It  is  expected  that  a  program  to 
produce  those  wing  parameters  described  by  Jones  will  fciiow 
from  this  development. 
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III.  aATRIX  NOTATION 


Though  egudtions  using  subscripted  variables  are  a 
uri'uiy  method  of  portraying  ma  theaiatical  relations,  it  may 
oj  >.!  :iei  to  V  Lsaali2e  what  is  happening  in  those  eguatiojjs 
^he  -  ']i  ?y  o  ri  written  in  matrix  format.  In  this  chapter 
egu  iM.  -  of  reference  1  which  lend  themselve  to  tliis 
type  of  notation  are  rewritten  in  matrix  format.  The 
subjcripts  (:ri3ices)  i,j,k  are  described  for  the  wing  in  the 
0e  I  lane  <;  s  rhaown  in  figure  3.  1. 

It  will  le  convenient  at  many  points  throughout  this 
presentation  to  transform  two  dimensional  matrices  into 
vectors  and  conversly  vectors  into  two  dimensional  matrices. 
For  doing  this  the  following  sytolic  relations  are  defined; 


Fig  i:  3.1  as  it  is  drawn  represents  a  two  dimensional 
matrix  witii  i Tices  (J,I)  and  the  index  K  shows  the  relation 
of  a  column  vector  that  has  been  transformed  into  a  matrix. 
The  inverse  operation  can  also  te  deduced  from  this  figure. 


A.  CIRCULATION  AND  PBESSUHE  EEIATIONS 

The  following  equations  are  the  general  circulation 
derivatives  and  pressure  relations.  Those  equations  taken 
from  LeftL*  i:c‘^  1  are  noted  by  tiieir  ejuation  number  preceded 
ly  Iitt*;  I.  Therefore  (J6.8)  becomes; 


Simrlurly  for  (J6.9) 


1  1 


"M*.  \ 


(3.2) 


[wj-([6j*[vj)[rj 


tvxl  tvW 

By  iefining  two  new  matrices  [A^  and  [B]  as  follows. 


(3.68) 


[a]  --  [  A']  +  [\/] 

t¥»6V 

[b] --[b']^[v] 

6V!(t¥  feVic^V  6V».6v 

the  final  wing  slope  functions  can  be  written  as 

[w,]^[a]  [r] 

tVxl 

for  the  additional  lift,  or  as 

6Vxi  4y/i4V  tVil 


(3. 69) 


(3.  70) 


(3.71) 


(3.72) 


for  the  basic  lift.  The  additional  lift  ejuation  can  bi 
invuited  as  follows  to  be  used  in  later  calculations 


(3.73) 


(M”''  CUl 
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[Q]^[H]fE] 


(3.64) 


s-<r 


Now  another  matrix  can  be  defined  which  shall  be  termed  the 
V  matrix.  natrix  V  is  a  partitioned  matrix,  which  is 
divided  into  three  parts  as  shown  in  figure  3.6. 


V  - 

_ 


Figure  3.6  Matrix  V. 


Therefore  the  final  form  of  the  integral  equation  is. 


(,fAl  feVxl 


3.  Fin^l  WiHi  Slc£e  Function 


(3.65) 


The  two  relations  for  Wp*  and  Wp”  can  now  be  super- 
imp  jred  to  jive  the  final  equation  for  the  wing  slope  func- 
tio:.  Fion  (.17.2) 


^  ^  < 
■  >  ■' 


•+  Wa 


(3. 66) 


After  sul  s  itu  t  i  ng  the  relations  for  Wp’  and  Wp"  and  simpi- 
fyi  .  ti  e  t  ]  i:itions  tecome 


'  * 

X  / 


fevliv  iyni 


(3.67) 


24 


(3.57) 


\ 

K I 


HtcI 

2.  Second  Winil  5 lo pe  Function 


The  integral  equation  for  the  second  wing  slope 
function  (J7.25)  is 


\yp  -  zTT/fi 


.^2 


j  ^  ci© 


(3.59) 


This  equation  requires  only  the  derivative  of  the  circula¬ 
tion  function  in  the  ©  direction  along  the  trailing  edge. 
Therefore  a  relation  involving  the  same  E  matrix  as  before 
can  be  utilized.  By  defining  a  function 


-  2Tr/fi  ^ 


(3.59) 


equation  3,58  can  be  written  in  index  notation  as 

8 


(3. 60) 


In  matrix  notation  it  becomes 


(3.61) 

^  y 

Equation  3.42  can  be  substituted  for  far  and  the  equation  now 
becomes 


(3.  62) 


[w."j  [rj 

bif-i-  Fx(r  yxi 

Then  substituting  equation  3.44  for  17  gives 

[vv^.  j  -  [h  ][£]  ^24  ^  IrdJJ 

S-At 

For  conveni  •  iico  an  additional  matrix  is  defined  as  fellows. 


(3. 63) 
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S(kpyk(?^l)E''(^0'hhQ) 

^  5  £  (l^0^2y  ho]  (3.  53) 

The  ''•ilcul.itLon  of  [A']  and  [B*  '  then  is 

/\  ( >: -  T(  kp,  kQ-s)  OaikQ-h,  ^<?) 

T {kp^hQ)D2(kQ,hQ) 

■^T  (kp,  ko^s)  ko) 

-  S(kp^ho-i)  12 "( ko -i,  ko) 

-SCkPj  ko)E''i^o,kQ) 

-S  (kp,  kQ*-i )  E*(kQ*hho) 

-SCkPko*2)  E*Lko*2,  ko)  (3.54) 

b'Ci^Kko}  -  T(kp  ka-j')  Dh(kQ  -8j  ko) 

T  (kp,  ko)  Dh(tio^ko) 

->■7  (kp  ko*8)  Ob(kQty^  ko  ) 

~S(  kp^ko-l)  EX^q-i,  ko) 

-S  (kp, ko)^ikQjkQ) 

-5(kp,kQ*l)  BXkQ*/,ko) 

-S(kP  ko*2)  EXkQ*2ykQ)  (3.55) 

The  various  factors  that  occur  in  the  above  equations  take 
on  definite  numerical  values  if  and  only  if  the  indices  KP 
and  KC  are  first  assigned  definite  numerical  values.  When 
this  is  done  the  above  equations  define  the  value  of  a 
single  element  of  either  [A']  cr  [  B*  ].  The  matrices  [A’] 
and  [B'  ]  can  be  stored  as  two  64  by  64  arrays.  So  the 
matrix  equations  for  the  first  ving  slope  function  are,  for 
the  additional  and  basic  lift  cases,  respectively, 

[v/;]-- fAj[rj  (3.56, 

6Ynl  6VK.I 
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the  integral  equation  can  be  written  in  index  notation  as 


The  atrix  D*  represents  either  Da**  or  Db*^  whichever  is 
appropriate.  The  elements  of  Da*  and  Db*  are  shown  in  table 
1  and  those  of  E *  in  table  2. 


To  make  the  distinction  between  additional  and  basic 
lift  the  following  twc  matrices  are  defined 

.'.V 

/V(  -Z  ( r(kf,k}D:(k,h)-S(hf,,k)E'(k,kQ})  (3.  »9) 

t-  / 
tv 

B  ( kpj  ^  ( r(  (k^kQ)-  5(kp,  k )  E  '‘(nM))  (3.50) 

The  above  caculation  for  [A']  and  [ B* J  is  simpler  than  the 
6U  by  64  multiplications  indicated.  If  one  uses  tables  1 
and  2  the  products  are  simplified  to, 

r  -  T(kp,kQ-S)  P*(Kci'^,kQ) 

+  T(  Kpj kQ) D*(kQ,  hkt) 

^  T ikp, kQ<-s}D^(kQ*f^  ka)  (3.51) 

T(kp,k)  Di(kj !?(?)  =  T  (kf/ kG>-S')  Ohikq-S'^  kcf) 

^Tikp,k(i)DUkq,kQ) 

^  Tikp,kq^y)0^\kQ^^,kQ)  (3.52) 

64 

E  S(kp,k)E.''(k^kiq)  -  Sikr,kQ-<)  ^"^{ko-l.kQ) 

^  E)ikpjha)E*ikq,ko) 


I 


Q,  .i)-A{^r.(ih2r,(7)'irr(s)] 


(3.41) 


Thi.'S  .':et  of  ejuations  can  be  written  as  a  single  matrix 
egu  it  ion 


A.  ^  c-  frvi 


8«5-  S*l 


(3.42) 


This  E  matrix  is  the  same  as  the  one  described  earlier  in 
figure  3.4.  From  equation  (J7.15)  at  the  trailing  edge 


Rewritten  as  a  vector  equation 


(3.43) 


rO’1))  rr(s7) 
r(5-6)J "  26  (  r(Q^) 

g*-\. 


(3.44) 


Notice  that  of  the  64  elements  in  the  P  vector  only  the  last 
16  are  used. 

C.  TNTECRAL  EQUATION  FOR  iING  SLOPE 
'ii.R'1  Slqge  Function 

T'lp  integral  equation  (J7.24)  ,  is  restated  here  as 
egua  t ion  3.45 


2  IT 


J  ^iNf  Fej  (3.45) 


By  defining  the  following  two  relations 


T(hp,k)- 


5//V© 


(3.46) 


20 


there  will  be 


on  the  values  ♦1/2A0and  -1/2A0.  Likewise, 

two  other  points,  at  locations  K=(KQ-1)  and  K=(K^*1),  at 
which  retakes  on  the  values  +1/24©  and  -1/2A0.  This  situ¬ 
ation  is  illustrated  by  the  schematic  diagram  in  figure  3.5. 


Figure  3.5  Circulation  point. 

More  complicated  situations  arise  if  the  point  Q 
lies  immediately  adjacent  to  a  boundary.  In  this  case  the 
number  of  field  points  having  ncnzero  values  of  and  f®  nay 
range  from  three  to  five.  The  details  can  be  seen  from 
careful  study  of  the  matrices  Da*,  Db*  and  E*. 

Along  the  trailing  edge  a  set  of  relations  for  the 
circulation  function  and  its  spanwise  derivative  is 
regiired.  From  a  similar  analysis  to  that  used  for  the 
spaiJise  l«:rLvative  over  the  wing,  the  following  set  of 
eguiMon  for  the  trailing  edge  is  obtained. 

\^.r(i)=  (3.39) 

r^r  (j)  =  ^[2l^(/'0*2rr(/*lj] 


1  9 


j=2,3,...,7  (3.40) 


In  a  manner  similar  to  the  equations,  equations 
3.31  through  3.33  can  be  writter.  in  single  index  notation 

re(k)--£[r(h)*ir(hO]  k=1,9,17 . 57  13.35) 

[^(hhA[ir(k-l)>ir(k<0]  k=2-7,10-15 . 58-63  (3.  36) 

. « 

These  equations  when  viewed  in  the  vector  form  are  much 
neater. 

The  'j-’iiGtaJ  c  .-'nfiguration  of  the  E*  matrix  along  with  the 
description  of  its  elements  is  shown  in  figure  3.8  and  table 
2.  Notice  from  the  figure  that  all  of  the  nonzero  eleinonts 
of  n.3trix  E  *  lie  along  just  four  principal  diagonals. 
Storage  can  be  minimized  by  storing  [Da*]  and  [Db*]  as  64  by 
3  matrices  and  [ E  *]  as  a  64  by  4  matrix,  capitalizing  on  the 
large  number  of  zeros  in  each. 

It  is  useful  to  consider  the  hypothetical  case  in 
which  the  circulation  function  is  assigned  unit  value  at  an 
arbitrary  point  Q  (the  'circulation  point')  with  matrix 
indices  J2  and  IQ  or  equivalent  vector  index  KQ,  and  is 
asst  ined  zero  values  at  all  other  points  of  the  calculation 
gr  i  1 . 

The  above  hypothetical  distribution  of  P  implies 
that  of  the  many  field  points  (denoted  by  index  K)  there 
will  be  certain  ones  adjacent  tc  the  circulation  point  Q  at 
which  the  derivatives  and  P*  take  on  nonzero  values.  In 
the  simplest  and  most  typical  case  where  Q  does  not  lie 
adjacent  to  any  of  the  boundaries,  there  will  be  two  field 
points,  at  locations  K=(KQ-8)  and  K=  {XQ+8)  ,  at  which  Intake:; 
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For  the  derivatives  in  the  spanwise(6)  direction  a 
set  of  relations  similar  to  those  for  can  be  shown.  From 
(J7.  17)  ,  (J7.18)  and  (J7.23)  respectively 


fl(  =  2^[r(  ^)]  (3.31) 

j  =  2,3,...,7  (3.32) 

12  a) -^6)  a)  -^rf7,l)-^Z  r(^y'^)l  (3.33) 


So  that  in  matrix  notation  the  equation  is, 

fc]^[E][r] 


and  the  E  matrix  is  defined  in  figure  3.4. 


f^{^,  yi)  -  ^  i“2,3,...,7  (3.27) 

r^(^,s)-MHr(f.7)^r3r(/,^)j  (3.28) 

for  which  the  matrix  notation  is; 


or  solving  for  circulation, 

S'nS'  g-xSf 

with  the  contents  of  Db  shown  in  figure  3.3 


(3.29) 


(3.  30) 


Figure  3.3  Ratriz  Db. 

Notice  that,  the  element  of  Db  are  identical  to  those  of  Da 
excii't  for  til*"’  first  two  elements  of  the  first  row,  which 
are  f  and  0  vice  1  and  1/3. 


16 


^J.e  <  t  Li;  ^jcntial  conf  iijurat  ion  and  description  of  the 
eieiiM  iits  II  t  .le  matrix  Da*  are  shown  in  figure  3.7  and 
table  1  rei'ptc  t  i  ve  ly .  Notice  from  the  figure  that  all  of 
the  nonzero  elements  of  the  matrix  Da*  lie  along  just  three 
principal  diagonals. 

Since  in  the  basic  lift  case  the  pressure  at  the 
leading  edge  is  finite,  an  additional  term  appears  in  the 
equation  for  the  circulation  derivatives.  That  term 
involves  the  second  derivative  cf  P.  Since  is  zero,  equa¬ 
tion  (Jo. 28)  for  the  pressure  at  the  leading  edge  is. 


i*!  8^  "S'  g-n  ( 


or  solving  forf^^^ 


(3.21) 


[  c.  j = f  r  nJ  [ft] 

6*1  fs*i  S'*/ 

Dy  taking  the  elements  of  both  and  {P,.} 

them  in  the  first  column  of  the  partitioned 
the  same  name  thusly. 


(3.  22) 

and  placing 
matrices  witli 


(3.23) 


[Kj=[{P]i[o)] 

yx/  '  gxT 


(3.24) 


the  vector  equation  can  then  be  written  in  matrix  form 

(3.25) 


S'xS'  Sxi  s^xfr 

Now  the  equations,  {J7.10),  (J7.  16)  and  (J7.14)  are; 


(3.26) 
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Note  that  this  eijuation  involves  the  transpose  of  both  the  T 
matrix  and  the  rj  matrix.  A  description  of  the  matrix  Da  may 
be  found  in  figure  3.2. 


Figure  3.2  Matrix  Da. 


Equations  3.13  through  3.15  can  also  be  written  in  a 
single  index  notation. 


1<k<8 

(3.17) 

9<k<56 

(3.  18) 

Qa)--rAr3r(k-s]^mk)l 

57<k<64 

(3. 19) 

In  order  to  get  adjacent  points  in  the  chordwise  direction, 
the  integer  8  must  be  added  to  or  subtracted  from  the  k 
index  describing  the  point  as  can  be  seen  in  figure  3.1  The 
above  equations  now  can  be  written  in  matrix  notation 

6VkI 


14 


At  the  leading  and  trailing  edges  (J6.15)  and  (J6.  16)  are 

ACpl}  -  Pl  (3.9) 

ACp)rJ  '  R-  (>)  -  I  Prj  (3.  ''0) 

^  g-Ki 

Generally  the  pressure  equation  (J6.22)  can  be  written; 

I'p  ] 


Srxjr  gxjt  gxS'  Sx? 

or  solving  for  the  circulation  derivative  in  the  chordwise 
(^)  direction 


*^3  K 


ir  1 

p  SIN^ 


sKi  S<i 


(3.  12) 


B.  CALCULATION  OF  DEFIVATIfES 
1 «  Chordwise  Derivatives 

Equations  (J7.9)  ,  (J7.  16)  and  {J7.  14)  are  written  in 
the  index  notation  as  follows 


[;(p)=r,[r(pi)^jr(p2)] 

(3.  13) 

[i  ^  ^ 

L  =  2  f  3  ,  m  m  ,  f  7  (3.14) 

(3. 15) 

This  helps  to  visualiae  the  equation  in  the  matrix  notation 
below 


[r,r-[o.ir7 


8^S'  g'xS 


(3. 16) 
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X  "  2  y 


h--<=i  h--l7 


h-z 


^'3  h  =  3  - 1 1 


M=^ 

ft  =  |2 

A  =  5 

wmm 

k-l^ 

IlL. 

k^l5 

ft;  5' 

k--l6 

ft  =57 


b--S‘=f 


b-60 


ft  =  6/ 


M-63 


k 


Figure  3.1  Index  notation. 


(3.3) 


Uouj  the  leading  and  trailing  edges  equations  (J6.  10) 
bhrvi  jh  {.J(  .  13)  respectively  are 

i  >^-'1  •  fr  ?  (3.“) 

S\l 

(  ii)  -  f  r  /  (3.5) 

l  I  Ur) 

g-x  I 

-  5  r  ?  (3.6) 

'  5X 1 

g-x/ 

The  I  ressuie  distribution  (J6.14)  is 


(3.8) 


2 


Elements  of  Da*  and  Db’* 
ilong  the  Three  Principal  Diagonals 


Numerical  coefficients  shown  must  be  divided  by  A© 


Table  2 
Elements  of  E*" 

Along  the  four  Principal  Diagonals  (cont*  d) 


Numerical  coefficients  shown  must  be  divided  by  2i© 


IV.  SUMMARY  OF  KEY  RELATIONSHIPS 

A.  PRELIfllNABY  CALCULATIONS 

The  initial  set  of  calculations  to  be  performed  by  the 
computer  are  for  the  values  of  those  constants  which  need 
only  be  calculated  once.  The  following  equations  use  only 
the  input  parameters  ,  X  ,  2r  and  the  indices  i  and  j. 

(4.1) 

(4.2) 

(4.3) 

(4.4) 

The  following  equations  involve  vectors  whose  indices  i  or  j 
go  from  1  to  8. 

(pU)^ 

cosCeij^)) 

yU  (^)  -  A"  ^ 

3  1 


(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
(4. 13) 


A0-  % 

A<^Ae 

C  -  2  TT/R 

^//  AQ 
C  =  Jrr^ 


(4.  11) 


0t.(/  ) 

0r(/'’ 

0T  ) 

©/vx(^  ) 


o 

7T 

O 


(4.  12) 
(4.  13) 
(4.  14) 


Indices  i  and  j  can  be  obtained  from  the  index  k  through  the 
two  following  integer  equations 


(4.  15) 
(4. 16) 


A.  =  (  ^/d)  I  (integer  division) 

^  -  k  ~S  ( 

Indices  ip  and  jp  can  be  obtained  from  kp  similarly. 

B.  CALCULATION  OF  THE  A  AND  B  MATRICES 
1  •  Al  and  B *  Matrices 

Since  the  matrices  [A]  and  [3]  are  the  sum  of  the 
[V]  and  [A']  or  [B'J  matrices,  each  of  the  matrices  will  be 
computed  separately  and  then  brought  together  to  form  [A  ] 
and  £B].  Recall  the  equations 


[,A']=[r][Q:)-[S)(£'] 

itf)=lrJ[D;]-fs]U-] 


(4.  17) 


(4. 18) 


Wher*-.  [T],  [0*]  and  [E*]  have  been  previously  defined. 

It  is  convenient  to  nake  the  following  variable 
definitions 
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(4.  19) 


-- 

/iX  =  X(A7,)*c5-  C/>C/> 

r,  ^VaX^^aV 

T (  kp,  k)  ~  C  (  “  Tp)siN  S(/) 

Sikp.k)'-  C'(^ 

T (kf,  kp)  -  C  ( ) SiN9(/p) 

S  (kr,  >!p)  -C'(  ^Y~t'  S,IN(t>iv) 


(4.20) 
(4.  21) 

(4.22) 

(4.23) 

(4. 24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 
(4.  29) 


A  cii^'ck  iihoul  1  be  made  to  see  whether  the  control  point  and 
iielu  point  coincide.  If  they  do,  then  equations  4.28  and 
4.29  sliall  l)r  used  instead  of  4.26  and  4.27  The  matrices 
[A']  and  [13']  from  equations  3.54  and  3.55  are 

f\(Mr.  ^<c?)  -  ^0-^)  0*(hQ-^y  ko) 

*T  (kryko)0*(ho^ko) 

*1  *sj  Out kq *ir, Po) 

-S(kP,kq-l)£-*(kQ-lykQ] 

-5ikP,kci)^*(hQyhq) 

-SCkpy  ko+l)  £  XkQ  tl^  kq  ) 

~S(kpJ^Q*2)E*(ko*2,ko) 


\ 
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B'(kP,h<Q)  -  Tikfi^kQ'S)  0^(kQ'?i,ka) 

{hf^  /?(?) OkikQ, kd  J 

*T(kP,  kQ*3)  Db(fiQ*Si  kn) 

-^ikp.kQ-i)  E-\kQ-i/ ko) 

-S(kp,kQ)E*(ku?,ka) 

-SC  kjo,kQ-fl)E*CkQH.kQ) 

~SCkP,kQ*2)  £'‘Ckq*2,kQ)  (‘♦.31) 

The  various  factors  that  occur  in  the  above  eijuations  take 
on  definite  numerical  values  if  and  only  if  the  indices  KP 
and  KQ  are  first  assigned  definite  numerical  values.  When 
this  is  done  the  above  equations  define  the  value  of  a 
sinjle  eJer'ont  of  either  [A']  cr  [B*].  The  matrices  [A*] 
and  i  B'  ]  c.'ii  be  stored  as  two  64  by  64  arrays. 

'  *  Ills:  1  Matrix 

Recall  figure  3.6  with  the  associated  equation  3.64 
rep.;-ated  ht  it;  to  provide  clarity. 

(4.32) 

(4.33) 

these  additional 


Z\Xr  '  \  id.yi,)-^c 

/  ^  ^ 
-V^Xr  Csri, 

AXr 

»r/ 


(4.34) 

(4.35) 

(4. 36) 

(4.37) 


[v]^[[o  ]  liWi'liW] 

(’‘it.i-i  iifY' 

[o]  --  [H]f  Ej 
fey  US'  fey.y 

In  addition  to  selected  definitions  above 
definitions  are  made 
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4 


(4. 38) 


(4.39) 

(4.40) 


A  cifck  j;hc'ul’.  be  made  to  see  whether  the  spanwise  coordi- 
iidt  .  coini.ido.  If  they  do  then  eijuation  4.40  shall  be  used 
instead  of  4.39  Then  [Q]  can  be  computed 

CH  h, 

Now  the  [V]  matrix  can  be  constructed  as  was  shown  in  figure 
3.6.  Since  the  first  48  columns  of  V  are  zeros,  only  the 
nonzero  elements  of  V  need  be  stored.  These  amount  to  a  64 
by  16  matrix. 

3*  1]}®  A  and  B  Matrices 

Finally  these  matrices  are  superimposed  to  give  the 
A  a:il  3  matrices. 


A(  V/  (kpyhoj-f-  V(kf^,ko)  *^2) 

kaj-B'C^p  ko)^V(kp,kQ)  (‘*•'♦3) 


This  gives  the  complete  A  and  B  matrices,  which  will  be 
utilized  in  the  following  calculations. 


C.  ADDITIONAL  LIFT 

For  Iht?  .specific  additional  lift,  as  explained  in  more 
detail  by  Joiils  [Ref.  1],  setting  the  normalized  slope  to 
-1,  gives 
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I  M 

(,  ■  '  I  1-1 

fcVxl 


(4. 44) 


'hi;:  circjl  ilinn  can  then  be  transformed  by  the  operation 
described  earlier  to  get 


Irj^ir] 


(4.45) 

J 

fevx)  S'*  Sr 

After  being  transposed  it  can  he  used  in  equation  3.  16  to 
get  the  chordwise  derivative 


r. 

■'xS' 


D^][f]  (4.  >16) 

S'kS'  S-xSr 

The  corresponding  derivative  in  the  &  direction  is  not 
generally  needed  and  is  not  shown  here,  but  has  been  stated 
earlier  as  equation  3.34.  Now  [ Q  ]^is  transposed  and  substi¬ 
tuted  into  equation  3.11  to  get  the  pressure  distribution. 


VxV  ^ 


(4.  47) 


lr>^y  yxir  yxfr 

Froni  this  pressure  distribution  one  can  calculate  the  aero- 
lynir  Lc  certtr  and  the  slope  of  the  wing  lift  curve. 


D.  l  ASIC  llfT 

Zn  tJio  case  of  basic  lift,  the  pressure  distribution  is 
initially  prescribed.  From  this  distribution  one  obtains; 


yxy  S-xir  s-xi-  SrxK 

and  additionally  along  the  leading  edge 


if 


(4. 48) 


(4. 49) 
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where  only  the  first  column  of  each  of  the  matrices  and 
[P^]  is  nonzero.  Ihen  and  [r,]  are  transposed  and 
substituted  into  equation  3.30  to  qet 


.r]'=[arrffr.,j‘fr.r 


(4.50) 


The  transpose  of  [fj^is  the  basic  circulation  which  after 
being  transformed  into  a  64  by  1  vector  is  utilized  in  the 
following  equation  to  get  the  wing  slope  function. 


(4.51) 


The  wing  slope  function  can  then  be  used  to  find  the  varia¬ 
tion  ,n  ta<'  l(?Lght  of  the  mean  wing  surface  above  the  refer- 
fc-nc>'  ':y  {  l.'.no. 
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V.  CONCIDSIONS 


From  the  eijuations  developed  and  listed  in  the  sumnary 
of  Key  relationships  of  the  previous  chapter  it  can  be  seen 
that  they  m<.y  be  translated  into  a  computer  program  w.hich 
cdljiiates  tJii'  wing  slope  function  and  the  lift  distribution 
of  1’/  dc  3:  J!  :i  ted  wing.  A  further  straightforward  extension 
of  Li-e  analysis  (not  given  in  this  thesis)  then  fixes  the 
various  aerodynamic  parameters  of  interest  such  as  the  slope 
of  tiie  wing  lift  curve,  the  location  of  the  aerodynamic 
center,  tlie  pitching  moment  coefficient  about  the  aerody¬ 
namic  center,  the  induced  drag  and  the  variation  in  the 
height  of  the  mean  wing  surface  above  or  below  the  reference 
xy  plane. 
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